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Abstract 

We extend our previous analysis of the toy model that mimics the mode coupling theory of 
supercooled liquids and glass transitions to the out of equilibrium dynamics. We derive a 
self-consistendt set of equations for correlation and response functions. 
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I. INTRODUCTION 

Recently we have introduced a mean field toy model that mimics the mode coupling the- 
ory(MCT) of supercooled liquids and glass transitions with trivial Hamiltonian [jl], ^, Analyses 
were limited to the equilibrium dynamics. An important feature of the model is that the strength 
of "hopping processes" 0] that destroys the non-ergodic state of the ideal MCT [||, |(| can be tuned 
so that nonergodic state is still allowed in some region of the model parameter space. This implies 
that the so-called hopping processes do not seem to be the same as thermally activated processes. 

In order to obtain further insights into the nature of MCT we consider the out of equilibrium 
dynamics of the model. In connection to this, recently the out of equilibrium dynamics of the mean- 
field-type spin glass models and other related glassy models was considered Q. In particular, the 
out of equilbrium dynamics of the spherical p-spin model §, represented by the closed set of 
equations for the off-equilibrium two time correlation function C(t,t w ) and the response function 
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G(t,t w ), was analytically solved in the long time regime P]. The analytic solution has revealed 
interesting features of out of equilibrium dynamics of the model. The system exhibits a strong 
waiting time dependence in the relaxaion of both C and G, i.e., aging behavior at low temperatures. 
Moreover, the fluctuation-dissipation theorem (FDT), i.e., the relationship bewteen C and G in 
equilbrium, is modified in an interesting way. Similar FDT violation have been observed in the off- 
equilibrium dynamics of supercooled liquids in computer simulations [1(J 11, 12 1 and an experiment 
II- 

We note that the all the out of equilbrium glassy features in the p-spin and related models are 
driven, as in the equilibrium dynamics, by the dissipative nonliearity in the equation of motion 
which comes from the nonlinear Hamiltonian. In the present toy model, as we see below, there is no 
dissipative nonlinearity since the Hamiltonian is trivial, i.e., gaussian without disorder. Instead the 
equation of motion involves the non-dissipative, i.e., reversible mode coupling nonliearities which 
drives the slowing down in the relaxation and the dynamic transition in the equilibrium dynamics. 
Our model possesses the reversible nonlinearities since we had a fluid in mind in constructing the 
model. Here we aim to see the out of equilibrium dynamics of the model driven by these reversible 
nonlinearities. 

As a first step in this direction we derive below the self-consistent closed set of equations for five 
correlation functions and five response functions. The method used is standard: the generating 
functional method in which two fictitious external fields are introduced for each dynamical variables 
entering the model []l4]]. But in view of more complications involving correlation and response 
functions of these variables we will sketch a derivation. We indicate possibility of reducing these 
complicated equations to the set of two correlation functions and two response functions. 



II. TOY MODEL 



Here we consider the toy model with M-component velocity-like b variables and iV-component 
density-like a variables, M being smaller than N. We have shown that in the limit of M, N — > oo 
with 5* = M/N finite, the parameter 5* becomes a measure of hopping Q. As a special case, if 
we take 5* = 0, then we obtain the zero-hopping model for the variables a, that is, the model is 
trivially non-ergodic. For 5* = 1 the hopping fully contributes and the system is always ergodic. 
For intermediate values of 5* we expect an ergodic to nonergodic transition at some value of 
T = T(S*). 

In this paper, after introducing the model we consider off-equilibrium dynamics that eventually 
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involves 5 correlation functions and 5 response functions. 

Our toy model is described by the variables a>j with j = 1, 2, • • • , N, and b a with a = 1, 2, • • ■ , M 
which are sometimes abbreviated as x which spans the phase space of the model. The model is 
described by the following Langevin equation 

LO 

hi — Ki a b a H j=Jij a CLjb a 
VN 



ba = -'yba - uj 2 K ja aj j=:J ija (uj 2 aiaj - T(%) + f a 

V N 



< fa(t)U(t') >= 2 1 T8 a p8{t - t') (1) 

where the f's is the thermal noise with zero mean and the angular bracket is the thermal average 
over such noise, and the usual summation convention for repeated indices are used. Here and after 
we will use Roman indices for the component of a and Greek for that of b. Here 7 gives a decay 
rate of the variable b a and u is seen to give a measure of the frequency of the oscillation of the 
variable a,-. 

For later purpose, we require that the matrix Ki a satisfies Ki a K^ = 5 a p. It is then easy 
to show that the equilibrium stationary phase space distribution of the Fokker-Planck equation 
corresponding to (Q) is given by [|], 

D e (x) = cst.e ^j=i 2T * Z - Q = 1 2T a (2) 

where est. is understood to be a suitably chosen constant. 

The mode-coupling coefficients JJy a are considered to be static random variables satisfying the 
following statistical properties: 



"ija — 0) 



2 _ 

(Sikfiji + hi5jk)5ap + Kip{K ka 5ji + K la 5j k ) 



g 



JijaJkl/3 — jy 



+ Kjp{K ka 5u + Ki a 5 ik ) 



(3) 



where 7TT ' 7 is the average over the independent Gaussian distribution of the J's. Eventually we 
take the mean field limit M,N — > 00, keeping 5* = M/N finite. 

In writing the model equation the temperature and other model parameters are fixed during 
time evolution, for example, in a situation after the quench. Naturally the model equation is valid 
only in such a time region. 
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III. ACTION INTEGRAL 



In order to analyze the toy model in the limit M, N — > oo, we introduce the following generating 



functional [14]: 



Z{h a ,h a ,h b ,h b } = J d{a} J d{b} J d{a} J d{b} 

x exp ji J dt(h a jaj + hfaj + h b a b a + h*S a )} 



x e So+Sl (4) 



where 



S = J dt \idi{di - K ia b a ) + ib a (b a + jb a + uj 2 K ia ai - (t) 

Sj = Jj ka Xj ka (5) 

Xjka = J dt ~^ {-idja k b a + ib a (oj 2 aja k - T5 jk )} (t) (6) 

In the above we have set the Jacobian of transformation of variables to unity assuming the ltd 
calculus M. In the limit M, N — > oo, we find that the last term —TSj k in the integrand of (||) can 
be neglected, and will be dropped from now on. That is, we take 

Xjka = J dt ~fjj {-ia-ja k b a + iiv 2 'b a aja k } (t) (7) 

We then notice that the replacements idj — ► (io 2 /T)aj, b a — ► b a /T on the rhs of (^) make this 
term vanish. Hence we can also rewrite (0) as 

Xjka = f dt ^j^ {-ia.ja k b a + iiv 2 b a aja k } (t) (8) 

where 

LO 2 1 

idj = idj + —aj, ib a = ib a + —b a - (9) 
The quantities of interest are the out of equilibrium correlation functions' 2 ] 

C a (t,t') = i < >, C ab (t,t') = ^K ia < ai (t)b a (t') >, 

Cfea(M') = -^r^ia < (0 >> CtftO = 7^ < b a (t)b a {t') > 

C* (tt') = -L < a* (i)a* (t>) >, a£ = K ja a 3 (10) 



[1] A consequence of choosing the ltd convention is the causality condition on the response functions. This implies 
that when the responses of a(t) or b(t) to the disturbances o(t') or b(t') occur simultaneously, the limit t' — ► t must 
be chosen in such a way that t is always greater than t' . 

[2] The definitions of C a b and Cta and those of the G's do not matter in the end. So we will use the symmetrically 
defined ones. 
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and the response functions 

G a (t,t') = — < aj{t)iaj{t') >, G ab (t,t') = ^K ia < ai (t)ib a {t') >, 

G ba (t,t') = ^K ia < b a (t)iai(t') >, G b (t,t') = ^< b a (t)ib a (t') >, 

Cf (tO = ^ < £ (*') >, ™* = K ja i aj (11) 

We note that the new types of correlation and response functions C^f and G^ are needed to obtain 
the closed set of equations for C's and G's for M < N. For M = N we have C^f = C a and 

U o — U a- 

Since we are here concerned with out-of-equilibrium situation we will not use the time translation 
invariance nor the fluctuation-dissipation theorem (FDT)^. We now take averages of (||) over f a 
and the J's where we use 

i J dtb a (t)f a (t)\ _ -jTfdtb a (t) 2 



Then we have 

<e S °>=e S °, e S i =e Sl (13) 



where 



50 = J dt{iai{ai- K ia b a ) + ib a (b a + jb a + u 2 K^ai + jTib a )^ (t) , 

51 = 9 - j ' dt dt'6(t-t% ja (t) 



iiija(t') + i jia (t')) 



(14) 



[3] The usual FDT takes the form where the response function is proportional to the time derivative of the correlation 
function. But here in the case of gaussian Hamiltonian the FDT is given by 



G a (t-t') 


= -6(t-t 


')^c a (t- 


-*'), 


G a b(t — t 


= -e(t-t')^c ab (t-t 


G ba (t-t') 


= -$(t-t 




-*'), 


G b (t-t') 


= -e(t-t')±c b (t-t') 


Gf(*-0 


= -0(t-t 




-t') 







where 8(t) is the usual step function equal to 1 for positive t and zero otherwise, the appearance of which comes 
from the causality. Another property arising from the causality plus the above FDT is the following for arbitrary 
X(t)=X(a(t),b(t),a(t),b(t))- 

< A(t)X(t') > = < X(t)A(t') >= for t > t' 

where A(t) — (a(t),b(t)) and the indices are suppressed for brevity. This fact is only limited to equilibrium when 
h = h = 0. Hence < X(t)A(t') >= for t > t' will not be used here. For h — and arbitrary h, however, the 
causality requires < A(t)X(t') >= for t > t' , which will be used later. 
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and 



£ija(t) = 4(-MffljM(t). CtfaC*) = ^(*6aOiOj)(t) 



(ib a ai(ij)(t ) 



(15) 



Here we remind that we can interchange and £i 3 - a in 5/, @, as we have noted in connection 
with (||) and (0). Therefore we can split 5/ into the 4 contributions 



5/ = 5f + Sf + + 5f . 



(16) 



where 



N 



sr = y - Ut dt'e(t-t')it ia (t) 



+ K if3 K ka (( kjf3 (t') + £"*./?(*')) + K jl3 K ka (^ kif3 (t') + iikp{t')) 



Sf = ^ I dt I -,')<;<,,,(,, 



A r 



s* = |- dt dt>e(t-t')£Ut) 



N 



sf>^ g - T jdtjdt'e{t-t')i\ ja {t) 



(§*«(*0+Sfa(*0) 



(17) 



Here the terms giving non-vanishing contributions in the equilibrium were overbraced. But this is 
no longer enough in out of equilibrium situation as we shall see. As an illustration we look at the 
first term of (|l7| ) or its integrand <Sj a . That is 

,2 



or 



-iai(i)aj(i)& Q (£)][-mj(tO%(OMO] 



w 2 A^5*iai(t)ia i (t / )C a (tt , )^(«') + 



(18) 



where the ellipsis contains quantities like < iaiit) • • • > which are absent for h, h = and also 
contains other fluctuation terms. Due to the presence of a factor A" in front of the second member 



of the rhs of (18), these fluctuation terms disappear in the limit M, N — > oo with a finite 5* . We 
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then analyze each factor like £,ij a (t)£,fj a (t') in ([l?]) in the limit M,N — > oo, where we only impose 
the causality condition t > t' reflected by 0(t — t') in ( |17] ) but not FDT. 

After tedious but straightforward algebra with h = we arrive at the effective quadratic action 
given below. To simplify the expression we introduce the following notation^: 



X Y(t) 



dt'X(tt')Y(t') 



In writing down the quadratic action below, for simplicity, we suppress time arguments and indices 
for the variables a, b, so that for instance we write K • b for Kj a b a in matrix notation. We also 
omit integral signs and tg) for the moment. The total action is given by the following matrix form 



Hot 



&eq ~t~ ^eq ~t~ ^oe 



( 



(fa, ih) ■ Cl 



eg 



ia | * /a 

+ (ia, ih) ■ Q eq 
ih I lb 



+ (ia, ih) ■ n r oe 



V 



(19) 



where [eq] and [oe] stand for equilibrium and off-equilibriun, respectively. Also note S oe is absent. 
We have 



eq 



n + 0/ 

/ 



On Onm 



LjV; 



TZ nh K 



(20) 
(21) 

(22) 



-66 



and furthermore, 



a 



eq 



iL eq 



n 



eq 



0° + Q J 

/ 
\ 



l N d t -K 
iv 2 K T l M (d t + 7 ) 

S6 a K T lM^bb 



(23) 
(24) 

(25) 



[4] Note signifies causality. Alternatively, this causality is taken care of by redefining the E's by absorbing 8(t — t') 
into them. Then the above integral is from — oo to oo. If desired one can do the same for the response functions. 
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and 



/ 



(26) 



liv(-S ao + AE ao - 2AE® ) -(S afe + AS afe )K 
^ (-E 6a + AS 6a - 2 AS® )K T -1 M (S 66 + AS 66 ) 
where 1n(1-m) are the unit matrix of rank N(M), K is the N x M matrix with the elements Kj a 
and K T its transposed M x N matrix. All these matrices are multiplied by a matrix whose tt' 
element is the delta function 5(t — t'). Also Onm, ®mn are the matrices of N x M, M x N, 
respectively. Here the memory kernel S's are defined by 



Sa6(tt0 = 



T 



(5*C a (tt')C b (tt') + (**) 2 C a6 (tt / )C 6 a(tt')) 



-2^-J*C a (tt / )C 6o (tt / ), S ba (tt') 



-2^<5*C a (ii')C a6 (tt') 



2<?V 



(27) 



The other types of kernels AS and AS® are defined as follows: 

AS aa = gW(5*G a C b + (S*) 2 C ab G ba ) 
AS® = g 2 uj\5*C a G b + {5*) 2 C ba G ab ) 
AS a6 = -g 2 u 2 6*[C a G ba + G a C ba ] 



(28) 
(29) 
(30) 
(31) 

2g 2 u 6 5*C a G ab (32) 
Note that every term of the AS's contains one factor of the G's and one factor of C's in contrast 



AS 6b = 2g 2 u>*C a G a , AS 6a = -2g 2 ujH*C ab G a 



AS® 



to the S's, (27), each of which contains the two factors of C's. 



IV. CORRELATION AND RESPONSE FUNCTIONS 



We now proceed to response and correlation functions. We first introduce correlation and 
response matrices C and G with sub-matrices C aa =< aa > and G aa =< aia > etc whose 
elements are < aa >u,jt'=< cn(t)aj(t') > etc. Thus the entire correlation and response matrices 
are written as 

(33) 



V 




Caa C ab 

C ba c bb 

The formal matrix equations determining correlation and response matrices take the form which 
are obtained from the effective action defined through (|H]) to (|26|) : 

n-G = 1 (34) 
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O • C = (£l eq + £ii q ) • Gt 



(35) 



(36) 



A. Response Function 



We first take up the response functions. The equations for them are written in terms of sub- 
matrices as follows where equilibrium and off-equilibrium parts are separated: 



(9t + E aa )G aa — (1 - E a fc)K • Gfe a 

+ [[(-£«, + AE aa - 2AE® )G«a - (E afe + AS ab )K • G 6a )]] 

= In (37) 

(dt + H aa )G a fe — (1 — S a f,)K • Gbb 
+ [[(-E oa + AE aa - 2AS®)G a6 - (E a6 + AE a6 )K • G 66 )]] 
= Onm (38) 

(cv 2 + E fca )K T • G aa + (8 t + 7 + E 66 )G 6a 
+ + AS fea - 2AE® )K T • G aa - (S feb + AE feb ) • G 6a ]] 

= Omn (39) 

(a; 2 + S fea )K T • G ab + (& t + 7 + S feb )G feb 
+ [[(-E ta + AE fca - 2AS® )K T • G ab - (E 66 + AE&&) • G w ]] 
= 1m (40) 



Here [[• • •]] are off-equilibrium parts which vanish in equilibrium due to the FDT which makes all 
the AE's and AE® equal to — E's. From this we can deduce the equations for 5 response functions 
in the following manner. 

We define the following notation valid for arbitrary matrix X, Y of ranks N, M, respectively: 



We first apply tr a • • • to @, next tr 6 K T • • • to (g|), then tr a K • • • to (f|), also tr b ■ ■ ■ to @, and 
finally tr b K T ••• K to @. 



We then end up with the following set of 5 equations for 5 response functions where again 
equilibrium and off-equilibrium parts are separated: 



^^E*^ tr fe Y 



M 



\(d % + E aa )G a - (1 - H ab )5*G ba ](tt') 
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+ [[(-E ao + A£ aa - 2A£® )G a - (S afe + AX ab )6*G ba ]}(tt') 

= 6(t - (41) 

[(ft + £ aa )G ab -(l-£ o6 )G 6 ](tt') 
+ [[(-S ao + A£ aa - 2A£ a )G ab - (S afc + AS a6 )G 6 ]](W / ) 
= (42) 

[(a; 2 + S ba )Gf + (d t + 7 + S 66 )G 6o ](W') 
+ [[(-E 6a + AS fca - 2AE® )Gf - (S fcb + A£ fcfe )G ba ]](tt') 
= (43) 

[(a; 2 + £ ba )G afe + (ft + 7 + S 66 )G 6 ](tt') 
+ [[(-E 6a + AS fea - 2 AS® )G afe - (S b6 + A£ 6b )G 6 ]](tt') 
= 8(t - t>) (44) 

[(ft + S oa )Gf -(l-S a6 )G 6a ](tO 
+ [[(-S ao + AS aa - 2A£® )Gf - (Soft + A£ ob )G 6a ]](iO 
= 6(t-t') (45) 



B. Correlation Functions 



We start with the formula (35) where the lhs are the same as those of fl37|) to (4C) except that 
the G's are replaced by the C's. Thus we need to consider only the rhs. Note that the rhs is the 
same as the equilibrium case if G^(ii') is given. G^(ti') in terms of submatrices is given by 

/ 

G\tt') = 



< ia(t')a(t) > < ia(t')h(t) > 

< ib(t')a(t) > < ib(t')h(t) > 



Also one can work out 



(Ci eq + £il g )(tt>) 



(46) 



j^[E«a(tt') + £«a(f*)]ljV, T[Z ab (tt') + ^ ba (t't)]K 

T[^Z ba (tt') + Z ab (t't)]K T , T[2j6(t - t') + E bb (tt') + Z bb (t't)]l M 

2^ aa (tt')l N , 2TS ab (tt')K 
^ 2TZ ab (t't)K T , 2T[ 1 5{t-t') + Y lbb {tt')}l M 

The last equality in (^7j) is due to the following symmetric properties of E's under exchange 
of two times: T, aa (t,t') = E aa (t',t), T, ab (t,t') = -4S 6a (t',t), and H bb (t,t') = T, bb (t,'t). These 



(47) 
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follow directly from C a (t,t) = C a (t',t), C b (t,t) = C b (t',t), and C ab (t,t') = C ba (t',t). The matrix 
equation (|35|) can be split into 4 submatrix equations as follows where the lhs is abbrevitated as 
d t C aa (tt') + ••• etc. 



d t C aa (tt') + ■■■ = 2— E aa (^) < ia(.)a(tO > 

+ 2T£ ab (t.)K- < »b(«)a(0 >, (48) 

d t C ab (tt') + ■■■ = 2-^E ao (t.) < i&(«)b(0 > 

+ 2TE a5 (i.)K- < ib(.)b(0 >, (49) 

dtC ba (tt') + ■■■ = 2TS ab (.t)K T - < ia(.)a(t') > 

+ 2T[ 7 <5(i - .) + S M (*.)] < ib(.)a(0 >, (50) 

d t C bb {tt') + ■■■ = 2TE ab (.t)K T - < ia(«)b(0 > 

+ 2T[ 7 <5(i - •) + E 66 (f«)] < ib(.)b(t') > (51) 



The rhs are the same with the equilibrium case. 

We are ready to find the rhs of the equations for correlation functions, which can be done 
following the same procedure as for the G's. Thus we first apply tr a • • • to (pj|), next tr b K T • • • to 
the rhs of (g|), then tr a K • • • to the rhs of @, also tr b ■ ■ ■ to the rhs of @) ; and finally tr ft K T K 
to the rhs of @. 

Results are summarized below in the form of 5 self-consistent equations for 5 correlation func- 
tions: 

l(d t + E aa )C a - (1 - Z ab )5*C ba ]{tt') 
+ [[(-E aa + AS aa - 2AE® )C a - (E afe + AZ ab )5*C ba }](tt') 
= 2?-Z aa (t.)G a (t'.) + 2TZ ab (t.)5*G ab (t'.), (52) 

[{d t + E aa )c ab -(i-E ab )c b }{tt') 



[[(-Eaa + AE aa - 2AE® )C ab - (E afc + AH ab )C b }}(tt>) 
2^ aa (t.)G ba (t>.) + 2TX ab (fr)G b (*'•), (53) 
(w 2 + E fea )Cf + (d t + 7 + E 66 )C 6o ](tt') 
[(-E 6a + AE fea - 2AE® )Cf - (S feb + . 

2TE a6 (.t)Gf (*'•) + 2T[ 7 <5(i - •) + E 66 (i.)]GU(^), (54) 



+ [[(-E 6a + AS fea - 2AE©)Cf - (E feb + AE b6 )C ba ]](tt' 
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[{uj 2 + ^ ba )c ab + (d t + 7 + z bb )c b ](tt') 

+ [[(-E 6a + A£ 6a - 2AE® )C a6 - (£ b6 + AS 6b )C 6 ]](tt') 

= 2T£ afe (.i)]G fea (t'.) + 2T[ 7 <5(t - .) + S 6b (^)]G b (t'«), (55) 

[(s t + s aa )af -(l-EaftjCfcjctO 

+ [[(-Eaa + A£ aa - 2AE® )Cf - (E a6 + AE afe )C fea ]](tt') 

= 2^E aa (i.)Gf(i'.) + 2TE^.)G a6 (i'.) (56) 



The 5 equations (|4l|)-(f45|) and 5 equations (|52"|)-(|5"6|) constitute 10 equations that self-consistently 
determine 5 correlation functions and 5 response functions. Note the rhs of d52])-(|5"6|) are the same 
as in equilibrium case. The lhs again have been split up into equilibrium and nonequilibrium 
portions as we have done in (|4l|)-(|45|). 

We now give several comments related to the set of equations (|52|)-(^6|): 



So far we have considered a general out of equilibrium situation where the model parameter 
remains constant. Now we specify the condition to meet a typical aging experiment. We then 
suppose that a system in some equilibrium state is quenched at the time t = and starts to 
evolve into another equilibrium state which characterizes the model parameters. We let the 
system age till some time t w (> 0), and we measure at a later time t(> t w ). This would mean 
that in the above equations we should take t' = t w and the time integrals denoted by • = s 
should be over the region s > 0, which is further limited by the causality conditions on the 
E's and the G's. Thus each integral in the equations of the response functions, (pl|)-(^5|), 
is in the interval t w < s < t. Similarly, each integral in the lhs of the equations for the 
correlation functions, (|52])-(|56[), is in the interval < s < t whereas each integral in the 



rhs is in the interval < s < t w . We then observe that the rhs of (52) to (56) are the 
source of contributions to the correlation functions from thermal noise generated after the 
quench. Note that the functions that multiply the Cr(t'»)'s in the rhs of (|52|) to ( |56| ) are the 



correlation functions of renormalized thermal noises. See the equations (25) of ||. 

In the discussion above we have not considered effects of the initial condition at the time of 
quench, say, to ]l6|]. This can be studied by inserting the properly normalized weight factor 
proportional to exp(— H/ksTo) where H is the system Hamiltonian and To is the initial 
temperature. For spin glass cases including Potts or p-spin systems, H contains quenched 
random interaction parameters, which have to be included in averaging the exponential of the 
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action integral over such quenched random parameters contained in the action integral. This 
results in additional terms in the final effective action. In our toy model, this complication is 
absent because the Hamiltonian is free from quenched disorder. The initial condition enters 
only through the initial values of the correlation functions, that is, C a (t = 0, t w = 0) = Tq/lo 2 , 
etc. 

• If we take all the AS's and AX®equal to — S's, consequences of FDT, we recover the equilib- 
rium equations. In this case the rhs of 5 equations (|5^)-(|56|) must vanish for t > t' because 
of the FDT which tells that the G"s to be proportional to the C"s. The direct verification of 
this follows if we note that in equilibrium t w can be shifted to 0. This makes the integration 
interval of s, and hence the rhs of (^2|) to (|56| ) to vanish. 

• No AE® fe appears in constrast to AE^, which is due to asymmetric way a and b variables 
enter dynamics. 

• No C^f and appear in the S's. One can see this by inspecting the derivations of (|l9|) to 
(|26|). We see that no terms containing the combinations 

a a( t ) a a{ t ')^ a a (*') w ith = Kj a aj, ia^ = Kj a icij appear. 



V. CONCLUDING REMARKS 



In the foregoing sections we have derived the exact self-consistent equations for correlation 
functions and response functions for our toy model which, in some sense, is complementary to the 
works of Latz [17] who studied the microscopic fluid system. 

The set of ten equations of our self-consistent scheme would be too complicated for further 
analyses. Now, the velocity-like b- variables can be made rapidly decaying by choosing a sufficiently 
large value for 7. In particular, this gives rise to a possibility of adiabatically eliminating the 
velocity-like variables as we have done to derive a Fokker-Planck type equation for the probability 
distribution function of a- variables only in |@, 
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